% Ejemplo de invocaci??n
% rk4(0,48,10,0.01,'f1');
% como ultimo parametro recibe la funci??n del modelo
% estan disponibles f1 y f2 (el que tiene u)

function energy=Puntod1()
% Parametros de Runge kutta-4
a = 0;
b = 48;
ya = 10;
h = 0.01;
func = 'f1';

% modelo de simulaci??n
kArr=[0.25,0.5,1,2,5,10,15,20];
global ku;
global k;
global TOUT_M;
global TD;
global Tstar;

TD=20;
TOUT_M=10;
ku=0.7;
colors = {'1', '2', '3', '4', '5', '1*', '2*', '3*'};
labels = {'k=1/4','k=1/2','k=1','k=2','k=5','k=10','k=15','k=20'};

% Curvas de la temperatura T(t) para distintos valores de k
for i=1:length(kArr)
	k=kArr(i);
	RK = rk4gral(a,b,ya,h,func);
	T=RK(:,2);
	Tstar=(-k/ku)*(TOUT_M-TD)+TD;
	energy(i,1)=k;
	energy(i,2)=Trapecio(a,b,h,RK,'Energy1');

%plot(RK(:,1),RK(:,2), colors(i))
%legend(labels(i));
%	hold on;
end

plot(energy(:,1), energy(:,2), '-r*');

%title('Energia para los distintos valores de k (primer modelo)')
xlabel('k');
ylabel('Energia');

end

% Rotulos del grafico
%title('Temperatura del Sistema T(t) para distintos valores de k y Temperatura Exterior Tout(t)')
%xlabel('Tiempo (s)');
%ylabel('Temperatura (?C)');
%axis([0,50,0,40]);

%plot(Tout(RK(:,1)),'2+');
%legend('Tout(t)');


%hold off;
